Insect pollinators can unlock an annual monetary value of more than US $100 million from crop production in Rwanda

Insect pollinators provide a natural ecosystem service to more than 80% of known flowering plants, many of which are part of our diet. However, their importance in Africa and an agriculture-dependent country like Rwanda has yet to receive attention. This encumbers policy formulation and investments in insect pollinators as a strategic agronomic input. Meanwhile, Rwanda cultivates crops that can benefit significantly from insect pollinators for superior agronomic outputs. To uncover this, we characterized the dependence of the crop production subsector on insect pollinators. Using the bioeconomic approach, we assessed the total economic value and the value due to insect pollinators of crops cultivated in Rwanda. We also evaluated the crop's production value per ton and whether production would meet consumption demands in the complete absence of insect pollinators. Using 71 representative crops currently grown in Rwanda, we found a direct dependency of 62% on insect pollinators. Of 32 representative crops used for economic valuation in two years (2014 and 2020), their total monetary value is estimated at $2.551 billion to $2.788 billion. Direct insect pollinator-dependent crops accounted for 20% (2014) to 18% (2020) of this value, with the share attributed to insect pollinators above $100 million. The sector's vulnerability to insect pollinators decreased from 7.3% in 2014 to 4.3% in 2020. The mean production value per ton of the direct insect pollinator-dependent crops was found to be higher in 2014 before declining in 2020. Using 21 representative crops from 2014 to 2020, we found that many direct insect pollinator-dependent crops will struggle to meet consumption demands in the complete absence of all suitable insect pollinators. Finally, we propose interventions and future research that could be undertaken. These insights are a critical first step to propel the government to act on insect pollination to support its food security agenda.

Insect pollination is one of the imperative ecosystem services for more than 80% of known flowering plant species, many of which are part of our human diet 1,2 .These pollinators differentially improve crop quality/quantity and the availability of viable superior seeds 3 .Studies show that 75% of leading insect pollinator-dependent crops could account for 35% of global food production 4 .When the world was disaggregated by its economic regions, 59.8% of crops in developing regions were found to depend directly on these pollinators 5 .From a financial perspective, at the global scale, 51% (n = 46) of crops with a direct insect pollinator dependency accounted for 39% (625 billion euros) of the world's production value (€1618 trillion) 6 .The share from insect pollinators was estimated at €153 billion, while the agriculture sector's vulnerability to insect pollinator losses was about 9.8% 6 .
Food production constraints and economic losses due to declines of suitable insect pollinators could increase our dependence on already extensively consumed carbohydrate-dense staples (cereals, tubers, pulses), many of which are wind pollinated 5,[7][8][9] .Insufficient consumption of a diverse number of nutrient-rich foods increases our vulnerability to non-communicable diseases [10][11][12][13] .Annual global deaths of up to 691,000 due to inadequate intake of insect-pollinator-dependent fruits, vegetables and nuts, which provide essential vitamins and oils, have been postulated from modelling studies 14 .Disease and death burdens would be higher in upper-middle-income countries 14 , while production losses would be higher in low-income countries, several of which are in Africa 14,15 .

Methodology
In this study, we assumed that all pollinators behave similarly (completely absent or completely present).This is irrespective of which crops they pollinate or how many pollinators can pollinate a given crop.Secondly, we used the mean pollinator dependency range/ratio of crops and not their actual pollinator dependency values.Data manipulation, analysis and visualization steps were performed using the R platform and suitable packages (Table_5_SuppInfo) 36 .

Data extraction
Crop price and agriculture variables (yield, production and area harvested) data for the period 2010-2020 were downloaded from the Food and Agricultural Organization website 37 in 2021 and then updated in February 2023, just before the completion of this study.These were downloaded as single files and then merged into one data set.Additional crops mentioned in the grey and published literature were used to enrich this data set.This included crops captured in the Agricultural Household Survey by the National Institute of Statistics of Rwanda 32 , the National Agricultural Export Development Board annual report 33 and the Baseline Report on Rwanda Horticulture 38 .We extracted distinct crops from this data set, which resulted in a total of 71 representative crops.This data set was used to map the insect pollinator dependency of these crops.Crop type information was extracted from the Indicative Crop Classification (I.C.C.) list 39 .

Crop-insect pollinator dependency mapping
The crop-insect pollinator dependency classification model by Klein et al. 4 was used to classify the above crops.The insect pollinator dependency depicts the extent to which a plant depends on insect pollination for improved agronomic outputs (yields) or the production quantity the plant would lose in the complete absence of suitable insect pollinators.These dependencies are represented by lower and upper limit values between 0 and 100%, referred to as the insect pollinator dependency ratio, rate or range of a crop.
The model classifies plants into three main groups: direct, indirect, and none insect pollinator-dependent. Based on the above values, direct insect pollinator-dependent crops are further divided into four subgroups: essential (> 90-100%), great (> 40-90%), modest (> 10-40%), and little (> 0-10%).Indirect insect pollinatordependent crops are divided into two subgroups: crops whose seed yields increase in the presence of insect pollinators and we consume their resultant vegetative parts (increase seed production, ISP), and crops whose seed production increases only during breeding, but we do not consume their resulting vegetative parts (increase breeding, IB).The latter crops reproduce through vegetative propagation.None insect pollinator-dependent crops do not require insect pollinators (wind-pollinated).Crops like tree tomatoes, stevia, moringa, and leeks with an unknown insect pollinator dependency were excluded.
The above classification model was further confirmed against country-specific crop dependency data sets from Brazil 40 , Argentina 41 , India 23 and Burkina Faso 29 .Where discrepancies occurred, we used information from the most recent publications or the African-specific context (Table_1_SuppInfo).Stenchly et al. 29 used a different classification for the indirect dependent group.They classified carrots, taro, onions and cabbages as directly dependent on insect pollinators, given that we consume their vegetative parts, whose production is influenced by seed quality/quantity achieved during insect pollination 42,43 .
The insect pollinator dependency table was generated after assigning an insect pollinator dependency (mean ratio, subgroup, group), crop name, species name, and crop type to each of the 71 crops in our data set (Table_1_ SuppInfo).Crops with a different dependency value defined by Stenchly et al. 29 were marked with an asterisk (*).

Evaluation of the crops' distribution
The percentage distribution of a given pollinator-dependent group, subgroup or crop type was computed using data in the crop-insect pollinator-dependency table (Table_1_SuppInfo).Firstly, we aggregated (subaggregated) crops based on their pollinator-dependency groups (direct, indirect, none), subgroups (essential, great, modest, little, IB, ISP, none) or crop types (fruits, vegetables, nuts/oilseeds, stimulants, pulses, tubers, cereals, sugars, spices/herbs).For a given number of crops in a subaggregate (n subaggregate ) and the total number of crops in the aggregate (n aggregate ), the percentage distribution of the subaggregate is computed as in Eq. (1).
For instance, to compute the percentage distribution of the direct insect pollinator-dependent group, we divide the number of direct insect pollinator-dependent crops (n subaggregate ) by the number of crops in the main aggregate/entire data set (n aggregate ) and multiply by 100.

Determination of the economic value of insect pollinators to crop production
The bioeconomic approach was applied to determine the economic/monetary value of insect pollinators to crop production 6 .It takes into account three input parameters: total crop production quantity (Q), farm gate price (P) and the insect pollinator dependency value (D).Hence, for a given crop (i), where i Є [1: I], the economic value (EV) is the product of the price (P i ) and production quantity (Q i ) as in Eq. (2).The economic value due to insect pollinators (IPEV) for a crop (i) is the product of its EV and the mean of its dependency ratio (D i ) as in Eq. (3).The final results are expressed in the currency of the price.
To calculate the EVs, we downloaded the FAO data set on the value of agricultural production from 2014 to 2020 as a single file.It contains crop data and the calculated economic value (U.S. dollars; $; Eq. 2) for several crops.We combined this with the data set on crop price, yield, production, and area harvested (Sect.2.1) to easily compare different variables and obtain additional crop names, price and production data.Where applicable, we used the price and production information to calculate the EV of crops as outlined in Eq. (2).Crops lacking price and production data, or an EV, were eventually excluded, resulting in a final data set of 32 crops.
Using data in our crop insect pollinator-dependency table (Table_1_SuppInfo), we assigned an insect pollinator dependency (mean ratio, group and subgroup) primarily from Klein et al. 4 as a new data set or from Stenchly et al. 29 , as a separate data set.We added the species name and crop type to these crops before calculating the IPEV (Eq.3).All downstream steps were applied to both data sets.However, for visualization, the detailed outputs from the data set of the former authors are often displayed and described in the results.
For visualization, we aggregated crops based on their insect pollinator dependency group and crop type before computing their total EV, IPEV, area harvested and production for each year.Then, we displayed the results for 2014 and 2020 (Table_2_SuppInfo).The annual outputs mean and standard deviation for the entire study period (2014-2020) were also displayed based on their insect pollinator dependency groups when the classification model of Klein et al. 4 (Table_3_SuppInfo) and Stenchly et al. 29 (Table_4_SuppInfo) was applied.

Determination of the vulnerability of the agricultural sector to insect pollinators
We calculated the vulnerability of the crop production subsector (CVR) if faced with a complete pollinator loss 6 by dividing the IPEV sum by the EV sum for any crop category.The result is multiplied by 100 and expressed as a percentage, as in Eq. ( 4).In our study, the IPEV for indirect and wind-pollinated crops is zero when D i = 0. Hence, this vulnerability solely depends on direct insect-pollinated crops when D i > 0.
We computed the vulnerability ratio for each crop type and the annual vulnerabilities using the total economic values generated in Sect.2.4 (Tables 2, 3, 4_SuppInfo).

Determination of overproduction based on consumption demands
To understand whether production can meet consumption demands without insect pollinators, we used the food balance sheet data from FAO for Rwanda for 2014-2020 6 , which contains consumption and production data.Crops and their products (e.g.maize and its products) were first aligned to the main crop (maize).Processed foods (e.g.infant food), crop types (e.g.cereals and others) or miscellaneous were excluded as they do not represent a single crop.This resulted in a final data set of 21 crops.These crops were assigned an insect pollinator dependency (mean ratio, group and subgroup), species name and crop type as defined in our crop-insect pollinator dependency table (Table_1_SuppInfo).
For overproduction before pollination loss for a specific crop (i), where i Є [1: I], its consumption (C i ) is subtracted from its total production quantity (Q i ) and then divided by its consumption as in Eq. ( 5).For overproduction after pollination loss, we obtained the proportion of production that is not dependent on insect pollinators (1-D i ), multiplied this proportion by its total production quantity (Qi), subtracted its consumption (C i ) and divided the result by its consumption as in Eq. ( 6).These outputs were expressed as a percentage of consumption.
Q i = total agriculture production (tons) of a given crop i, C i (tons) = consumption of crop i, D i = the pollinator dependence rate of crop i (%), (1 − D i ) = portion of crop i without insect pollination.

Determination of relative changes
Relative changes between two years of interest were calculated using the method described by Azien et al. 5 and information in Table 2_SuppInfo.For this, one of the years was set as a reference year.Therefore, to compute the relative change in a variable (x) between two years, the value of a variable (x) in the reference year (x t reference ) is subtracted from that of the year of interest (x t ) before dividing the output by the value of variable (x) of the reference year (x t reference ) as in Eq. ( 7).In our study, 2014 was set as the reference year.Our variables included yield, production, area harvested, production value, economic value and consumption demand before or after the complete loss of insect pollination for any given crop.The final outputs are expressed as percentages.
With x ṫ = value of a variable (x) for a given crop in year t, x treference = value of variable (x) for a given crop in the reference year.

Missing data imputation
The Multivariate Imputation by Chained Equations (MICE) package in R was used to fill in missing data 44 .MICE has the advantage over other methods of filling in missing values many times based on the observed values of a given sample and the relationship in the data set 45 .The MICE imputation was used to fill the production value of green beans (2014 and 2015) and green chillies/peppers (2014).

Statistical analysis
Statistical differences in production value between the three insect pollinator-dependent groups were calculated using the non-parametric Kruskal-Wallis test for each year.This is recommended for non-normally distributed data, especially if the researcher chooses not to manipulate the data through log transformation or removing outliers, for instance 46 .This test was selected over the analysis of variance (ANOVA) since the data for specific years (2014, 2015, 2016, 2020) arguably satisfied the test for normality, even though they satisfied the homogeneity of variance assumption.Secondly, not all crop types are represented in each insect pollinator-dependent group-this narrowed the number of independent variables that could be included.
The Kruskal-Wallis test was performed using its built-in function in RStudio.We defined two hypotheses, null and alternative, and then used the probability (p) value to statistically interpret the results 47 .For the null hypothesis, H 0 implies that the means of the dependent variable (production value) is not significantly different

Dependence of edible crops on insect pollinators
Of the edible crops grown in Rwanda, 62% directly depend on insect pollinators We found that 62% (n = 44) of representative edible crops (n = 71) that are cultivated in Rwanda are directly dependent on insect pollinators (Fig. 1, table).Of these, the majority (14 crops) have a modest dependency, 12 crops a little dependency, 11 crops a great dependency and seven crops an essential dependency (Table_1_Sup-pInfo).Crops that did not require pollinators constituted the second largest group (21.1%, n = 15), followed by 16.9% (12 crops) with an indirect dependency.Among the latter group, the majority (n = 9) are pollinated by insects during crop breeding, whereas for the others (n = 4), insect pollination is vital for increased seed production.

Fruits, vegetables, and oilseeds substantially contribute to the direct insect pollinator dependency
The most representative crops were fruits, vegetables, nuts/oilseeds, and spices/herbs.The first three contributed the most to the overall dependence of the direct insect pollinator group.All represented nuts/oils have a direct dependency, while up to 77% (n = 17) of fruits and 58.3% (n = 7) of vegetables also have a direct dependency (Fig. 1, white bars).However, most nuts/oilseeds (75%) are at the lower end of the dependency grid of insect pollinators (modest -little), in contrast to 58.8% (n = 10) of fruits that occupy the higher end of the dependency grid (great-essential).For vegetables, there is a comparable distribution among the insect pollinator-dependent groups.
Among the least represented crop types, a majority have a direct dependency.We found that direct insect pollination is vital for 75% (n = 6) of spices/herbs.Of the six pulses, 83.4% occupy the lower end of the insect pollinator dependency grid.Of the two stimulants, there was an equal share, with coffee having a modest dependency and tea having no dependency.Tubers populated the indirect insect pollinator group, predominantly for breeding purposes.All representative sugar crops and cereals are wind-pollinated (Table_1_SuppInfo).

Economic value of crop production and insect pollination
The total economic value of crop production increased in Rwanda from $2.551 billion (2014) to 2.788 billion (2020) The total economic value of crop production from the 32 representative edible crops in Rwanda was in 2014 estimated at $ 2.551 billion (Table 1), with a production of 6.6 million tons across an estimated area of 2.1 million ha (Table_2_SuppInfo).Of this, the direct insect pollinator-dependent group had the highest number of representative crops (43.8%).This accounted for 20.4% ($521.5 million) of the economic value.Furthermore, its production (14.1% of the total) was slightly lower than that of the none-dependent group (15.7%).This former group also occupied the second-largest harvested area (31.8%).
The total economic value of crop production in 2020 increased by 9.3%, equating to $2.788 billion (Table 1).This was the second-highest economic value attained between 2014 and 2020 from representative crops (Table_3_ SuppInfo).There was also a 33.1% increase in relative production (8.8 million tons) and a 5.8% increase in harvested area (2.2 million ha), despite the pressures imposed by COVID-19.Additionally, the economic value of the direct insect pollinator-dependent group decreased slightly by 3.4%, resulting in an estimated value of 503.5 million.Nevertheless, this was still higher than that of the none-dependent group ($414.6 million; relative increase = 7.3%) but lower than that of the indirect dependent group, which still made the highest contribution of 58.5% ($1.869 billion; relative increase = 13.8%).
The direct insect pollinator group had a 4.5% increase in production (966,356 tons) and a 25% increase in its area harvested (826,475 ha) in 2020.Its production remained slightly less than that of the non-insect pollinatordependent group (1 million tons) and far less than that of the indirect-dependent group (Table_2_SuppInfo).This latter group accounted for up to 77.5% (6.8 million tons) of the overall production, even though its area harvested had declined by 15% (811,063 ha).
At the crop type level, indirect insect pollinator-dependent fruits and tubers contributed the most to the total economic value (Table 1).The respective amounts were $879.1 million and $737.2 million in 2014, before decreasing slightly in fruits (0.8%) but increasing in 2020 by 27.8% in tubers.Cereals, pulses (direct) and vegetables (direct) followed suit, with 130-354 million.

The economic value of crops due to insect pollination declined by 35.8% in 2020
The share of the total economic value of production due to insect pollination was in 2014 at an estimated $186.2 million (Table 1).This translated into a 7.3% vulnerability of the crop production subsector to insect pollination.

Table 1.
Total economic value of production (EV) in Rwanda in the year 2014 and 2020, share from insect pollination (IPEV), percentage relative change (RC) in the year 2020 against the reference year 2014, and vulnerability ratio (%VR).Source: Own calculations using FAO_data.M = million, K = thousand, USD = United States dollars.Insect pollinator dependency values (D i ) used to compute the IPEV were obtained primarily from Klein et al. 5 1), with a reduced vulnerability to insect pollinators of 4.3%.Stimulants (76.6%), fruits (71.9%) and vegetables (38.3%), represented mainly by coffee, mangoes, and pumpkins, contributed to this decline (Figure_2_SuppInfo).Amongst the crop types, vegetables (direct) accounted for the highest economic value due to insect pollination in both study years (Table 1).This value was in 2014, $153 million, before declining by 38.3% ($94.5 million) in 2020.They were also the most vulnerable crop type in 2014 (78.2%) and 2020 (68.7%).The pulses had an estimated value of 12.2-15.2million.Unlike vegetables, pulses were the least vulnerable to insect pollinators (5%).Fruits, estimated at $10.8 million in 2014, were the second most vulnerable crop type (38.7%).They also experienced one of the highest declines (71.9%) in 2020.The stimulant, coffee, was the fourth highest contributor in 2014 ($5.4 million) but experienced a 76.6% decline in 2020, resulting in nuts/oilseeds taking this spot.Stimulants also had the same vulnerability (25%) as spices, which made the least contribution in both years (Table 1).
When we recalculated the IPEV using the classification model of Stenchly et al. 29 , the total monetary value was estimated at $219.91 million in 2014 and increased to $235.24 in 2020.The sector's vulnerability was estimated at 8.62% in 2014 and 8.44% in 2020 (Table 4_SuppInfo).

The insect pollinator-dependent group had the lowest production per hectare of harvested area but the highest production value per ton
We also observed that in terms of production quantity per hectare of harvested area, the direct insect pollinatordependent crops were producing, on average, the lowest quantity per hectare (977,737.69tons/733,687.5ha = 1 .33tons/ha) over the 7-year (2014-2020) period (Table 3_SuppInfo).The wind-pollinated group had a slightly higher production per hectare of harvested area (934,986.5tons /514,047.43ha = 1.82 tons/ha), while the indirect dependency group generated the highest quantity (5,509,515.06 tons/875,269.86ha = 6.29 tons/ha).
The direct insect pollinator-dependent group had the highest mean production value per ton in 2014 and 2020.In 2014, this value was 769.9 ± 91.3$/ton (± standard error) and was significantly greater (p = 0.025) than that of the indirect dependent group (366.2 ± 51.5$/ton) but not the none-dependent group (584.0 ± 207.5 $/ton; p = 0.18).This value decreased in 2020 by 13.0% and 11.5% in the direct and indirect insect pollinatordependent groups, respectively, but increased slightly (5.2%) in the none-dependent group (Fig. 2).The value in 2020 remained significantly different (p = 0.010) only between the direct and indirect insect pollinator-dependent groups.Dry peas (none) had the highest production value/ton in both years, followed by groundnuts, coffee, lemon, oranges, mangoes (direct dependencies) and rice (none) in at least one of these years.

Insect pollinator-dependent crops could face major constraints in meeting consumption demands
Using data from 2014 to 2020, we found that among the eight representative direct insect pollinator-dependent crops, only coffee, and in selected study years, chillies, groundnuts and dry beans could meet consumption demands (Figure_2_SuppInfo).In groundnuts, this was apparent during the later years of the study (2017-2020), while in chillies, this was between 2015 and 2019.In contrast, oranges, soybeans, lemon/limes, and tomatoes had differential deficits ranging from 30 to 70% during the study period.The latter was affected after the loss of insect pollinators, whereas the others were affected under both situations.On average, over the seven years, oranges and soybeans showed deficits during both pollination situations, while the others had marked deficits (tomatoes and oranges) or declines (chillies and dry beans) only after pollinator loss (Table 2).Of the eight indirect insect pollinator-dependent crops, only potatoes, onions, and cassava could not meet consumption demands (Figure_3_SuppInfo and Table 2).While this was marginal in onions and cassava (± 10%), it was pronounced in potatoes during the early study years (2014-2018), with deficits of at least 40%, before recovering to exceed consumption demands (10%) during the later years.
Wheat and rice were the two wind-pollinated crops that could not meet consumption demands (Table 2), with deficits of at least 30% (Figure_4_SuppInfo).Although dry peas initially had a slight deficit (~ − 10%), they recovered in the subsequent years, remaining at marginal levels (~ 10-15%).Maize and sorghum, on the other hand, remained in excess of consumption demands throughout the study period, with excesses of about 30-100% observed.

Discussion
Our study has provided for the first time a more advanced insight into the dependence of the crop production subsector in Rwanda on insect pollinators, the economic value of crop production and the impact on consumption demands before and after the complete loss of suitable insect pollinators.What is unique in our study is the multiyear approach applied to observe fluctuations, if any, using the most recent data.
In our study, we found support for our hypothesis that insect pollinators are important to Rwanda's agriculture and will bring in a valuable monetary contribution.This was despite the absence of the added economic value of 30 insect-pollinated crops (e.g.passion fruits, macadamia, sunflowers, etc.) with no price data 6 .Losing this monetary contribution has repercussions on farmers' revenue 11 and GDP 48 .
Even though food choices are not guided by pollination, farming households in Rwanda (80%) rely heavily on their farm produce for consumption 10,32 .This includes insect-pollinated produce, especially vitamin A-rich vegetables (pumpkin, carrot, squash, amaranth/dodo), fruits (passion fruit, avocado, papaya), but also nuts/ oilseed (groundnuts), pulses (dried beans), spices/herbs and coffee 34,50 .Additionally, their consumption, which increases with availability, can be enhanced by promoting kitchen gardens and capacity-building programs on food and nutrition 13,51 .These actions align with the government of Rwanda's strategic plan to sustain a nutritionsensitive agriculture sector which increases dietary diversity and the consumption of nutrient-dense foods 34,35 .A focused intervention should, therefore, be on insect-pollinated vegetables and fruits, which, in our study, were the most vulnerable to insect pollinator losses and the highest contributors to the vulnerability of the production subsector.The latter vulnerability has remained slightly lower than that for world agriculture 4 but is higher in given years than for East Africa 6 .Together with nuts/oilseeds (soybean and groundnuts), these vegetables (tomatoes) and fruits (oranges, lemon) are among the crops that struggled to meet consumption demands in the 7-year study period.
Low-vulnerability crops (nuts/oilseeds, pulses, coffee) should also be given attention due to their added nutritional value (oils, proteins, vitamins) and international traction (coffee) 14 .To close nutrition gaps, biofortification is used to augment the nutrient composition of consumer-preferred crops (e.g.dry beans for iron delivery) in Rwanda 9 and globally 54 .However, if insect pollination is vital to these biofortified crops, appropriate interventions should be prescribed to achieve superior yields.For instance, dry beans can further benefit from yield increases per hectare of up to 46% during open pollination by insects 55 .
Our findings are consistent with postulations that low-income countries are yet to reap the full benefits of insect pollination 5,14,15 .Our representative direct insect pollinator-dependent group generated the least production per hectare of area harvested.Given that agricultural land is a limited resource 37 , Rwanda might struggle to compete in production quantity.However, they can capitalize on improving quality through strategic insect pollinator interventions as an added agronomic input.This will better position their produce in the international market.

Conclusion and recommendations
Given the potential value-addition of crop-insect pollination to food systems, nutrition and health, Rwanda should incorporate relevant insect pollination actions in its strategic plans.
A first recommendation is for Rwanda to increase efforts on beehive establishment, mainly where pollinatordependent crops are cultivated.This will directly benefit several vital crops like coffee, sunflower, macadamia, mango, avocado, etc., that depend on bees and other pollinators 4 .
Secondly, policies guiding bee hive establishment should align with those regulating pollutive fertilizers, herbicides and pesticides 56 .While these are strategic farm inputs in Rwanda 35 , their negative consequences on bees and other pollinators are a concern 20 .Insect-friendly farm inputs (organic fertilizers) should be promoted, while pesticides should be carefully analyzed for their impacts on insects and discontinued if harmful.
Thirdly, Rwanda should nurture natural ecosystems, which are pollinator-rich habitats 57,58 .Insect-pollinated crops could be cultivated close to such vicinities to facilitate interactions.Additionally, Rwanda's land consolidation action 35 should promote the establishment of patches of fields with high pollinator attractants (e.g.sunflower or flowering plants with bright colours) close to areas with pollinator-dependent crops 58,59 .
Above all, advocacy and capacity-building activities should be intensified 35 .Additional emphasis should be on the adequate handling of approved farm inputs, techniques to mitigate the impacts of harmful farm practices, the importance of pollinators in food systems, monitoring and documenting insect visitors, hand pollination techniques, and climate change mitigation actions (e.g.afforestation) 60 .

Challenges and future perspectives
This study has also identified gaps, mainly around the availability of data.Research efforts should be geared towards characterizing the insect pollinator dependency of crop cultivars in Rwanda, their nutritional value, and possible health benefits.This will open the gateway for impact-driven crop pollinator studies.A second research effort should focus on monitoring insect pollinator species in different ecosystems to clarify concerns on whether, how and at what rate insect pollinators are declining.Thirdly, national surveys and data collection efforts should capture crop production, area cultivated, area harvested, yield, farm gate price, and consumption demand data more consistently for an increasing number of crops, including minority crops in Rwanda.This data should also be disaggregated per crop and not aggregated under a crop type.These findings and recommendations could contribute to Rwanda's poverty alleviation and food security agenda, and could be translated in other African countries.

Figure 1 .
Figure 1.Dependence of 71 representative edible crops cultivated in Rwanda on insect pollinators.The cropinsect pollinator dependency ratio was inferred mainly from Klein et al.4 and other country-specific publications (Table_1_SuppInfo).The table in the chart depicts crops in each insect pollinator dependency group (direct, indirect, and none).The bar graph and associated digits represent the percentage distribution of each crop type.

Figure 2 .
Figure 2. Production value per ton ($/ton) of representative crops cultivated in Rwanda in 2014 and 2020.The values are represented with respect to their insect pollinator-dependent group (direct, indirect, and none).Grey dots represent the mean production value per ton in each group.Kruskal-Wallis and Dunn tests were used to calculate significant differences between groups using their functions in RStudio.
p > 0.05) among the independent variables (pollinator-dependent groups).Otherwise, we accept the alternative hypothesis H 1 if at least one group mean differs significantly (p < 0.05) from the others.Lastly, the dunn_test function in RStudio was used to obtain information about which two combinations of the insect pollinatordependent subgroups were statistically different.

consumption) Before. pollinator.loss After. pollinator.loss
This continues to be important in combating the 35% food insecurity risk, 33% stunting, 45% mortality risk (especially among children under five years), 26% obesity in women, and inadequate consumption of diverse crop types in Rwanda34,52,53 .